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Abstract 

How stable vortex-like magnetization in magnetic nanodisks with small aspect 
ratio {L/R « 1) is affected by two antidots is investigated analytically and by 
Monte Carlo simulations. For suitable ranges of the physical parameters this 
vortex presents bistable states when pinned around the antidots. The hysteresis 
loop obtained shows a central loop associated to the pinning mechanism of the 
vortex at the antidots. Our results agree qualitatively well with those provided 
by experiments and micromagnetic simulations. 



Recently, ferromagnetic disk-shaped nanostructures with submicrometer lateral di- 
mension (thickness) L < l/zm have been fabricated and investigated for their potential 
applications in a number of magnetoelectronic mechanisms. In particular, it has been 
observed that above the so-called single-domain limit, magnetic vortex states appear 
in these samples, exhibiting a planar-like arrangement of spins outside the core, where 
a perpendicular magnetization is observed[lJ. As long as one could manipulate these 
states other possibilities would emerge. In fact, one way towards this control is obtained 
by removing some small portions of the magnetic nanodisk, in such a way that the cav- 
ities (antidots) so created work by attracting and eventually pinning the vortex around 
themselves [21 [3l HJ [5l [6l [7j . Based upon such an idea, Rahm and coworkers [8] have 
studied the cases of two, three and four antidots (each of them with diameter ~ 85 nm) 
inserted in a disk with diameter ~ 500 nm, separated by around 150 nm — 200 nm. 
Their experimental results confirmed the previous statement about vortex pinning and 
put forward the possibility of using these stable states as serious candidates for mag- 
netic memory and logical applications as long as we could control vortex position, 
for example, applying a suitable external magnetic field which should shift the vortex 
core from one defect to another, and vice-versa. Basic logical operations have been 
obtained by means of bistable magnetic switching[^. Although experimental results 
are provided for this system a suitable analytical analysis is still lacking. The latter is 
important, for instance, to provide the basic physics behind the mechanism of vortex 
pinning and switching processes, giving the relevant parameters for a better control of 
possible applications and also indicating the limitations. Our present analytical model 
and Monte Carlo (MC) simulations have been able of capturing the basics involved in 
this problem, namely, how the vortex experiences the effects of the defects and how 
hysteresis loops are sensitive to the latter. 

Let us starting by considering a magnetic dot represented by a small cylinder of 
radius R and thickness L (so that its aspect ratio L/R << 1). In addition, we shall 
assume that along the axial direction (2;-axis), the magnetization M is uniform. Fur- 
thermore, if we introduce N isolated holes (each of them with height L and radius 
p « R) in the dot, the total magnetic energy of the nanodisk can be approximated, 
in the continuum limit, by 
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where A is the exchange coupling, D is the area of the cylinder face, m = M/Mg 
is an unity vector describing magnetization along D (with Mg being the saturation 
magnetization), = hmijn) = Hm/Mg is the demagnetizing field, /igxt is an applied 
magnetic field (Zeeman term) and /i = 1,2. The potential U , in turn, brings about 
the effect of the antidots distributed throughout the nanodisk, say, Hili Ui{r — di) = 
Ui{f- di) U2{r- d2) ■ ■ ■ UN{r- d^), with 

U,{f-d,) = { : ■ (2) 
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Therefore, the system of a dot with isolated antidots may be viewed as a cyhnder of 
radius R and thickness L with smaller cylindrical cavities with radius p << R, each 
of them centralized at di. Here, we shall study explicitly the case N = 2 (the treatment 
for N > 2 may be performed in the same way). Experimental and numerical results 
are available for very similar systems, say, disks with L ~ 30 nm, R ~ 500 nm whose 
antidots (up to four) have diameters 2p ~ 80 nm separated by 150 nm — 200 nm [S]. 

Now, let us consider a cylindrically symmetric vortex-like magnetization throughout 
a dot, say, with the vortex core centralized at r = 0. This is the ground state of a 
nanodisk in the absence of holes and external magnetic fields. For that, it is convenient 
to write rh = (sin 6* cosy?, sin 6* sin v?, cos 6*), with 9 = O^ir) and ip = arctan(t//x) ±tt/2. 
The function 6^{r) may be approximated by sin 6y{r) = in the dot center (r = 0), 
while sin 6'„(r) 1 far away the center, \ f \= r >> a {a = ^/A/Ms is the unit- 
cell element size or the exchange length; for most magnetic materials a = 5 — 6nm). 
In words, the magnetization consists of a small core where spins display out-of-plane 
components for regularizing the exchange energy, and an outer region where spins 
are practically confined to the dot plane face. In this case, the magnetic superficial 
charges in the lateral face of the dot and the magnetic volumetric charges (V ■ m) 
identically vanish yielding no contribution to hm- The antidots affect this picture as 
follows: from the point of view of the exchange term, antidots lead to a less exchange 
energy for the vortex. Then, the topological structure is attracted by the hole suffering 
modifications in its profile. In addition, the distribution of magnetic charges throughout 
the internal edges of the cavities (holes) and mainly along the external lateral face of the 
cylinder (whenever the antidot is not at the geometrical center of the dot) increase the 
magnetostatic energy due to a change in the product m ■ -hs {fig are unit vectors normal 
to external lateral surface of the disk and internal surfaces of the cavities). Thus, the 
demagnetizing field, hm, can be obtained from its associated potential = $v + 
"^edge + ^edge' usual Way, = — V$m- Here, "I'v is the magnetostatic potential 

related to the volumetric charges , while $edge $edge comes about from the surface 
charges on the external and internal (holes) edges, respectively. The contributions of 
the volumetric potential can be neglected since the approximations considered above 
leads to V ■ m = 0. 

For simplicity, the antidots centers are assumed to lie along a straight line that 
crosses the origin of the dot, say, at di and d2, as shown in Figured! We also assume 
that the vortex displacement / from its equilibrium position is not too large, so that 
it experiences no appreciable change in its profile ('rigid' vortex behavior ^lOj). Thus, 
the exchange potential experienced by the vortex may be estimated as [11] 

V^s, di, 4) = ^ In [(1 - s') (1 - /i - /2 + fj,)] , (3) 

where s = l/R, supposed to be small (| s \« 1), measures the relative shift of the 
vortex from the dot center; the functions fi are defined as = a^/(| di — sR p +&^), 
where b = 1.147p is a constant introduced to avoid spurious divergences whenever the 
vortex is centralized at one of the defects O [7]. Note that each fi is related to the 
attractive potential that each isolated hole induces on the vortex, while the product 
/i /2 accounts for a competition between them Thus, if our system were large 



enough, the vortex would displace towards one of the holes as shown in Ref.[TT]. In 
a small magnet, things are much more interesting once the vortex experiences a mod- 
ification in its profile and magnetization is no longer cylindrically symmetric. Then, 
the magnetostatic energy increases and a restoring force appears in order to pull the 
vortex back to the dot center. Indeed, for small displacements of the vortex we may 
estimate this energy shift analytically, like below 0, [ID] 



in which the first term is related to the contribution of the disk envelop, while the 
remaining ones are associated to the pinning of the vortex by one of the defects. In 
addition, Fi(^) = ^^^^ ^1 — ^—fj—j dt, where Ji{t) is the Bessel function. The 

possible three distinct situations are: i) the vortex is centered at the antidot 1, then 
a = 1 and /3 = 7 = 0; ii) it is centered at the antidot 2, thus a = 7 = while /3 = 1; 
iii) the vortex is not pinned at any defect, what yields a = /5 = and 7 = 1. 

Figure [2] shows how the effective potential, V^g = + Knag, behaves as function of 
s. Bistable states corresponding to the vortex-antidot pinned configurations appear for 
suitable values of the parameters. Our analytical model predicts that three parameters 
are the most relevant in this case: the radius of each antidot p, their center-to-center 
separation, Z) = — ^2!, and the (characteristic) exchange length a = \/~A/Ms. We 
have observed that the appearance of bistable states in V^fj are generally associated to 
larger values for p and a and to smaller as illustrated in Fig. [21 Actually, it shows 
that our analytical results qualitatively agree with those obtained in experiments of 
Ref . [8] , but there the bistable states are observed even for larger separations (~ 200 nm) 
than those fitted by the present analysis (~ 130 nm). This discrepancy is related to 
the small displacements of the vortex core assumed by our model. 

In the presence of the two holes the vortex equilibrium position, sq = {xo,yo), can 
be easily determined by evaluating dVes/ds = 0. Clearly, depending on the relevant 
parameters, the geometry of our system implies that sq will always correspond to 
(I Xo |> 0, 7/0 = 0) so that the general expression for the local magnetization along the 
dot face now reads: 

, . TV , . ±(a: - SqR) 

mx[x,y,Sn) = — , my[x,y, Sn) = — , (5) 

where the upper (down) signs are associated to counterclockwise (clockwise) magneti- 
zations. Since the system is antisymmetric under refiection against x-axis, the average 
magnetization along this axis vanishes, while the y component may be estimated like 
below (the factor 2 in the 2nd term accounts for 2 antidots): 
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In addition, if an external homogeneous magnetic field is applied along the ?/-axis, 
/lext = H /\Ms\ = hyj), then, in the lowest order, Ves must be augmented by : 

V^^Js) = -nM^hy (R' - p') (s - so) , (7) 

where s is the shift in the vortex center position caused by the field. So the total 
potential now reads V^otai = + H^xt • the presence of /lext the vortex equilibrium 
position is now Sh ^ sq. For example, if one takes the counterclockwise magnetization 
(upper signs in eq. ([5])), then for hy > {< 0) the vortex will be shifted to the left 
(right) of Sq. Thus, the equilibrium position of the vortex center in the presence of the 
applied external field, Sh, is simply a small shift of Sq along the x-axis. Therefore, in 
this situation (M^) = while {My) is calculated from eq. Q) with my{x, y, Sh) instead 
of my{x,y, Sq). How magnetization behaves as hext is varied gives the hysteresis loop 
left by the vortex motion under such conditions. 

Figure [3] displays the hysteresis loops for some values of the relevant parameters 
when an alternating-like magnetic field is applied along the y-axis. The central loop 
shown in this figure, on the left, is in qualitative agreement with the experimental 
results of Refs. [8], [9] . In these works, such a mechanism was observed for larger cavities 
radii and separations (p ~ 40 nm and D ~ 200 nm). In our analysis, we could observe a 
similar fact only for much smaller values of these parameters. Again, such a discrepancy 
comes about once our model is strictly valid for small s. The central loop, taking 
place at magnetic fields weaker than that for vortex annihilation, manifests the most 
important physical feature of the system: the existence of two metastable vortex states 
with the equilibrium position of the vortex center pinned at each antidot. There, the 
jumps are due to vortex core switching from one state to another while the plateaux 
corresponds to pinned vortex states. The mechanism of switching can be explained 
in the following way: in the situation shown in Fig. [2t^a) the vortex core can be 
initially pinned in equilibrium at one of the antidots, for instance, at the left antidot. 
However, this configuration may be perturbed if a strong enough external magnetic field 
is applied, say, along the y-axis, so that the vortex core is shifted to the another antidot. 
Inverting the magnetic field the vortex comes back to the initial configuration, but 
leaving a hysteresis loop in the complete round. The plateaux observed in the central 
loop of experimental data has a pronounced inclination respective to the magnetization 
axis, which may be associated to the deformation of the vortex profile on the dot 
face. Consequently, our preceding approach could not fit this fact. At the attempt of 
understanding this mechanism, we have performed MC simulations in order to improve 
the range of validity of our former results. For that, we consider the xy-model on a 
finite flat disc (of radius R) supplemented by a strong enough surface border anisotropy. 
Although it cannot be considered as the actual magnetostatic energy we expect that 
it could imitate its role, say, increasing the total energy as long as the spins develop 
normal components at the borders of the nanodisk and of the holes (2nd and 3rd terms 
below; see Ref. [12])). Indeed, this anisotropy has the effect of keeping the original 
cylindric-like profile of the vortex near these borders although allowing it to deform in 
other regions along the dot face. Therefore, Hamiltonian reads: 



H = -jJ2iS:S^ + SfS]) + B Yl iSk.rhkf + B Yl (^h.m,)\ (8) 

i,j kadiskborder hGholesborder 



where J > is the exchange ferromagnetic integral, Si = (Sf , Sf, Sf) are the classical 
spin vectors specified at the lattice sites i and the first summation is over nearest- 
neighbor spins. The other term represents the surface anisotropy where the sum over k 
includes only the border sites of the disk and sum over h includes only the border sites 
of the holes. Consequently, the unit vectors rhk are perpendicular to the circumference 
envelop of the disk and the unit vectors rhh are perpendicular to the circumference 
border of the holes. There, i? is a constant of single-ion surface border anisotropy so 
that for i? > (< 0) the border spins Sk and Sh tend to lie perpendicular (parallel) to 
rhj. and rhh direction respectively. Of course, here we chose the case -B > 0. Among 
other features, such a term enables vortex deformation on the dot face whereas keep 
its configuration at the border (for a fixed -B), so that we can now fit the inclination 
of the central loop discussed above. 

In our MC simulations for hysteresis loop we have adopted Metropolis algorithm [13] 
with the initial configuration of centered vortex in disks with diameters 2R = 40ao (ao is 
the distance between two spins in a discrete lattice inside the disk) . In addition, we have 
introduced two circular holes (regions with vacancy of spins) where we chosen radius 
p = 2aQ and center-to-center separation D = IOoq, in such a way that p/D 0.2 is 
close to that from experiments (see Refs.[8l|9]). Besides, we have taken the temperature 
T/ J = 0.05 and the parameter B/J = 0.03. This is the critical value for B/J so that, 
above it, non-centered vortices are energetically favorable in these disks (see Ref. |12]). 
Figure H] shows a typical central hysteresis loop as the vortex center is switched from 
one antidot to another and back again (note particularly that its inclination is in good 
qualitative agreement with experiments of Ref. [8]). 

In conclusion, our analytical as well as MC calculations have shown how two holes 
incorporated into the body of a magnetic nanodisk attract the remanent vortex-like 
magnetization to their centers, creating the possibility of bistable states of vortex- 
hole pinned configurations, as observed in experiments. The analytical model for the 
effective potential (exchange -|- magnetostatic) have yielded to this picture in good 
agreement with experimental findings. However, concerning the hysteresis loops, the 
analytical calculations show only a qualitative concordance with experiments. Indeed, 
the appearance of a central loop in those curves, related to the switching of the vortex 
core between the two stable states is verified here for values of p and D very smaller than 
those considered in experiments. In other words, our results in this case would lead 
to a good agreement with experiments only for very small displacements of the vortex 
core around its equilibrium position, Sh ('rigid' vortex regime). It implies in small 
distance D between the holes in a nanodisk. For improving this scenario (considering 
as much as possible generic distances) we have also performed Monte Carlo simulations 
which take into account the deformation of the vortex and effects due to the border. 
Our simulations for the hysteresis loop are in agreement to the experimental results 
and exhibits bistable magnetic switching mechanisms. The approach developed here 



also predicts the possibility of gyration of the vortex core about one of the three 
equilibrium positions (see Fig. [2]) with characteristic frequencies that depend on the 
relevant parameters R, p and D. The theoretical study of a quite recent observed 
gyrotropic frequency of magnetic vortex around antidots in nanosized structures'^ is 
under investigation and will be communicated elsewhere. 
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1 Figure Captions 



Figure 1: Top (left) and lateral (right) views of a nanodisk with radius -R, 
containing two antidots lying along the x-axis, whose centers are di and d2 apart 
from the disk center. 

Figure 2: Typical plots of V^s/A as a function of s. Here, we have taken 
R = 250 nm, L = 30 nm, p = 43 nm, and a = ^/A/Ms = 5.7nm (values considered 
in Ref. [9]). The defects are placed along x-axis at the positions: (solid curve) 
di = d2 = —50 nm; (dashed curve) di = —d2 = —80 nm. Note that the bistable 
states tend to disappear as long as the distance between the cavities becomes 
larger. 

Figure 3: The central loop associated to the switching of the vortex between the 
two antidots (bistable states) is jeopardized as long as D increases and/or a gets 
lower. External loops (not depicted) were associated to the creation/annihilattion 
of the vortex, (a) \di — d2\ = D = 40 nm, p = 10 nm, and a = ^/A/Mg = 17nm; 
(b) D = 40 nm, p = 10 nm, and a = 5.7 nm. 

Figure 4: Hysteresis loop for a disk with two antidots with radius p = 2ao and 
center-to-center separation D = lOao obtained by Monte Carlo simulations. 
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